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Abstract 

We investigate the critical scaling of the spin-1/2 antiferromagnet on the square lattice in the 
easy-plane (XXZ) regime, via numerical measurements of the entanglement entropy constructed 
from the zeroes of a polynomial ring. We relate these results to conformal field theory predictions 
for the area law scaling of entanglement entropyin the vicinity of quantum critical points in gapless 
conformally invariant 2d quantum systems, and in gapless systems with finite Fermi surfaces. Our 
measurements are focused in three low-temperature regions at fixed XXZ anisotropics of A = 1.01, 
A = 1.78 and A = 2.0 which probe the quantum regime between the Neel point and Ising phase. 



1 



I. INTRODUCTION 



The spin-1/2 Heisenberg antiferromagnet on the square lattice has had a long association 
with high-temperature superconductivity [l| 2], but there has been almost an equally long- 
standing controversy surrounding the topological nature of certain of its ground states [sl-jGl. 
This question has become more pressing recently due to the discovery, via STM and neutron 
scattering data, that spin impurities in the underdoped cuprates apparently induce a direct 



second order transition between the Neel ordered state and a condensed Bose phase 



A 



9|(l0|. In 



fuelling further interest in mechanisms of the quantum criticality of edge states 
contrast, the behaviour of the one dimensional analogue of this system is largely settled. 
Many quantum spin chain models can be solved exactly, and there is a well-known direct 
mapping between chains at criticality and the SU{2) WZW model which leads to 
a distinction between gapped and gapless groundstates for integer and half-integer spin 
chains. What makes the 2d system inherently more involved is that whilst the Hilbert space 
dimension is an extensive quantity the Hopf invariant is also zero 12|], which leads to the 
emergence of two competing spin wave theories for the ground state (or more properly a 
chiral wave theory coupled via a Chern-Simons term in the continuum limit 13|). Since 
the hard-core bosons of the Jordan- Wigner approach will then necessarily carry fermionic 



statistics in the 



ong-distance (IR) regime of this system (in addition to the short- distance 



(UV) regime [141 14|) this has motivated various ideas of Bose condensation via skyrmions 

nnn 

for the ground state of this system in the quantum regime |5l] |6l]|l5l]. 
In this article we consider the anistropic lattice model defined via. 



where i and j are spin sites indices, J is the nearest-neighbour spin interaction coupling, 
A is the XXZ spin anisotropy parameter, and S' is a spin-1/2 operator represented by usual 
the Pauli matrices with spin components r = {x,y,z). At A = 1 the zero temperature 
ground state is Neel ordered [l|. By considering the convergence of the support of spin 

nn 

correlators as a function of L it has been identified in [16|[17| that there is also long-range 
order (LRO) for anisotropy parameter values of < A < 0.13, and A > 1.78 (in the easy- 
axis and easy-plane limits). At A ~ a second order transition is expected as a function 
of temperature, but this is then expected to become first order when quantum fluctuations 



suppress this LRO [18| , which has been confirmed by numerics in [19j [20| . Attempts to treat 
quantum fluctuations in the spin-1/2 Heisenberg antiferromagnet on the square lattice in 
the vicinity of A = 1 have focused on treating the Neel order parameter as the basis of a 



2l| 



221 a 



low-energy effective fleld theory, via the (2 + l)-dimensional a model {^[sl. In 
fugacity expansion on the dual is used to identify that a dangerously irrelevant monopole 
contribution will lead to a second order (quantum) transition from the Neel to the Valence 

n 

Bond Solid phases 

The focus of this article is universality of the transitions in the quantum regime of the 
spin-1/2 Heisenberg antiferromagnet on the square lattice. However, two factors make the 
direct numerical veriflcation of the quantum transitions in 2l| 2^ difficult. Firstly, it is 



unclear what the relevant order parameters should be [10| , secondly, although lattices serve 



as effective short-range regulators, lattice IR cutoffs arise from statistica 



makes the lattice Lorentz symmetry (of gauge theories) approximate 231]. Although there 



averaging, which 



are no explicit gauge fields involved the construction of the spin-1/2 Heisenberg antiferro- 
magnet on the square lattice in (1), because the Hopf invariant is zero and the Hilbert space 



dimension is extensive 



121], it is important that an exact U{1) Wick rotation exists between 



Euclidean and Imaginary time at the boundary of the system. This symmetry is described 
as an "emergent" gauge field when the system is both Lorentz and scale invariant 2l| 22]. 
Therefore, because the Lorentz symmetry is approximate on the lattice, it is more difficult 
to construct numerics which have the correct Hopf invariant. For example, it is not possible 
to guarantee that a global minimum can be reached in quantum spi n sy stems using DMRG 
via the standard approach of minimising the following expectation 24| . 



{mm ■■ G (c^)^^ (2) 

Quantifying the scaling behaviour of entanglement entropy is one means of directly es- 
tablishing the universality of quantum transitions (at least for quantum chains). This is 
because scale invariance also implies conformal invariance in (1 + 1) dimensions. The quan- 
tum transition of spin-1/2 chains is characterised universally by the central charge c, where 
c = 1 for free bosons and c = 1/2 for free fermions. This has been verified by a number 
of numerical studies [2^ j^^] j28| (although similar problems arise with ensuring the Lorentz 



invariance of numerical lattice methods 
Von Neumann entropy {29!, 



241]). This universal scaling is quantified via the 



5=^^log(0+O(e-^/«) (3) 

where the Von Neumann entropy 5 describes the bipartite entanglement between two 
subsystems, c and c are the holomorphic and antiholomorphic central charges of the corre- 
sponding conformal field theory at the boundary between the subsystems, ^ is the correlation 
length (of the finite temperature system), a is the lattice UV cutoff, and x is the Euler char- 
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acteristic which categorises the topology of the boundary between subsystems j30| 31| 
Unfortunately, there is no generalisation of the above c-theorem for quantum chains 
to higher dimensions, since the Witt algebra associated with the primary (chiral) fields (of 
which the central extension is the Virasoro algebra) does not map directly onto the confor- 
mal anomaly in higher dimensions. The above area law scaling is therefore not universal in 
higher dimensions. On dimensional grounds, the coefficient in the area law scaling for the 
spin-1/2 Heisenberg antiferromagnet on the square lattice goes as a^~'^ and is not universal. 



but t 



lere shou^ 



d also be a universal term area law for the 2d system term proportional to 

Since we are faced with approximate Lorentz symmetries in standard numerical lattice 
approaches {23!, and the non-universal area law scaling of the entanglement entropy of 



the spin-1/2 Heisenberg antiferromagnet on the square lattice 25|] 30|] 3J], we take a new 



approach to identifying the quantum transitions of this system (which we have also recently 



applied to the quantum spin chain in 35|). For this, we construct an exact polynomial 
expansion for the partition function from Wick rotating the transfer matrix elements of a 
standard nonperturbative (Quantum Monte Carlo) ensemble. This enables us to construct 
two exact nonorthogonal polynomial representations of the primary chiral field separated 
via a nonperturbative Chern-Simons term (which corresponds to non-analytic portion of 
support of this polynomial) by identifying the quotient of this polynomial expansion. 
Taking the infinite volume limit of the support of this formalism corresponds to a (^-function 
renormalization of the partition function, as we show, and the (central) extension of the 
polynomial can be determined from this limit. In this article, we generate standard ensembles 
at three points in the quantum regime of the anisotropic spin-1/2 Heisenberg antiferromagnet 
on the square lattice (A = 1.01, A = 1.78 and A = 2 at /? = 100), and determine the scaling 
of the Von Neumann entropy (identified from the density of the polynomial zeroes) associated 
with the mapping on to its (1 + 1) primary chiral field at criticality. 



II. PARTITION FUNCTION ZEROES OF QUANTUM CRITICAL POINTS 



In order to understand the homology of 
on the square lattice that is presented in 
function in the following Wick-rotated form 



he treatment of the spin-1/2 antiferromagnet 
I it is instructive to write the lattice partition 
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where An, Bn and Vn are defined as the compact and noncompact portions of the spin 
operators in (1) determined through the nonperturbative matrix elements A(n) and s is 
the Euclidean-time lattice site index. To define the indices s' and s" a state is taken from 
G (C^)*^-^^ and the order of these tensor products is interchanged, which is valid on the 
system boundary 12]. Locally, however, this creates a mismatch between Oi and 0j {i = 1,2), 
hence polynomial representations of An and will not be orthogonal, in general, only in 
the special case where 6i = (pi and the system is both Lorentz and scale invariant. 

A meaningful statistical ensemble for spin-1/2 antiferromagnet on the square lattice can 
be obtained numerically using the continuous-time Quantum Monte Carlo method [s^] by 
generating a Markov chain from importance sampling using the following local transfer 
matrices |38| . 



L,L,T 



P{i,j},{i+hj} P{i,j},{i,j+1} 



^= n 

i=l,j = l,t=l \ P{i+l,j},{i+l,j + l} P{i,j + l},{i+l,j + l} 



(8) 
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J, 



where ^ and ^ for 

Jz/J-xy = A, and Ar is the Euchdean-time lattice spacing. Discrete steps in Euchdean-time 
are exchanged in the continuous-time method with discrete spin flips, hence the Markov chain 
;hat is generated is ergodic in probability space (Imaginary time) but not in Euclidean-time 



39l |. The nonpertubative matrix elements realised in this Markov chain can be identified 



at each step t, and up to a normalising factor oi (3J (the inverse temperature and bipartite 
lattice interaction coupling) these matrix elements are defined via, 



(*) 











1 - Ki,, 



(11) 



/i 



1 — ^t,i,j 

\0 1 / 

The Wick-rotation of these elements is then defined by compactifying the Euclidean-time 
boundary conditions of the lattice partition function from (0,/5] ^ 6 = [—11,71], via the 
following definition of the trace of the transfer matrix P (defined as the transfer matrix for 
the whole of the volume). 



TrP = |](1 - Xk^){l - {AXkf), Z= fve P 
k=i 



(12) 



where k is the lattice site index on Gi ® 02, and A is the spin anisotropy parameter. This 
formalism relates the matrix elements in (11) to those in (4). The assumption for this Wick 
rotation prescription is that each local plaquette on z ® j can be smoothly deformed into 
via the spin anisotropy parameter A. This is not true, since any two states in the system 
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are not orthogonal in general. Therefore, we recover an exact expansion in (3 J from these 
operators that is only analytic up to (pi ® 02 (rather than 6i (g) 62). However, this allows us 
to probe the criticality of the primary chiral field of the nonperturbatively realised system 
that corresponds to this sing ularity. 

In the replica method 30|, the Von Neumann entanglement entropy of an n-fold Riemann 
sheet (such as for n = L"^) is defined by identifying the analytic continuation properties 
of the lattice partition function with respect to Euclidean-time, 

S=-^^^l(7^) (13) 



n-»i dn \ (Z) 

where Zn is the lattice partition function for a subsystem consisting of n disjoint unions 
of Z. We go slightly further than the usual replica argument for defining the entanglement 
entropy, in our approach, by constructing a polynomial ring for C^. This removes the 
assumption that the partition function has to be factorisable in the above denominator 
{4^. This is particularly important for the spin-1/2 antiferromagnet on the square lattice, 
because general states in this system are defined by two non-orthogonal polynomials, hence, 
the entanglement entropy defined in the above limit is degenerate (even though the branch 
points are resolved in R^). The relevant quantity for determing the entanglement entropy 
in our formalism is the logarithm of detP. An analogous expression to the replica limit is 
then obtained by analytically continuing the following function from large s to s = 0, 

EVlnA. (14) 

k=l 

where = (1 — Afc^)(l — (AA^,-)^). This amounts to the renormalization of the singularities 
of the support of the free energy of the partition function via a (^-function prescription, where. 



k=i ^ ' -^^ 



(15) 



However, although the logarithm of det P (free energy) is completely defined by the first 
derivative of this ^-function at s = 0, the corresponding free energy minima is not unique 
since any number of the higher moments of C,p{s) can also be nonzero. Hence, to uniquely 
determine the entanglement entropy a further prescription must be given for resolving the 
choices of branch in k for all of these higher moments at a different point on the fundamental 
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strip, namely at s = 1, 



= Ve k{9) \n^{e), A; = ... 
=1 Jo 



This prescription is implemented by constructing a polynomial in 6 consisting of k terms 
- the zeroes of this polynomial ring therefore correspond to elements of the quotient 
(rather than M^). Hence, the meromorphic convergence of the zeroes of this polynomial 
uniquely specifies the resolution of the branches in k, via, 

IndetP = - / A(^) In(^) dO S = - A(^)lnA(^) (17) 

Jo Jo 

which is the familiar Von Neumann form of the entanglement entropy. 

Practically, the support of the logarithm of detP can be investigated via the scaling of 
the zeroes density of the eigenvalue problem det(P — {(3JY) = 0, since the divergences of 
the logarithm of detP occur where this polynomial ring has its zeroes. The zeroes are ob- 
tained numerically by rootfinding the characteristic polynomial equation for this eigenvalue 
problem, where the characteristic polynomial coefficients are obtained from powers of TrP 
by using Newton's identities (35|||4ll|, 



L2 k~l 

Z = ^{ck){(3Jf^, fccfc + TrP^' + ^c„TrP^"" = 0, n = l...L'^ (18) 

k=0 n=l 

This polynomial zeroes approach is similar in spirit to Lee and Yang's determination of 



the zeroes of the Ising partition function 4J], although the singularities in our formulation 



correspond to elements of the quotient (rather than Z2). In both cases the polynomial 
zeroes are constrained to lie on the unit circle in the complex plane of the polynomial 
expansion parameter X. This is because the quotient of the polynomial representations 
of each model and the Z2 symmetry of their spin operators form a polynomial ring, C = 
]R[X]/X^ + 1. For the Ising model this constraint is an identity, and all the zeroes lie 
on the unit circle, but for our formulation this constraint only applies at criticality, since 
PGL2{C) ^ PSU2 = 5*0(3). Motivated by this, we apply a simple difference relation to 
define the zeroes density A(/5J) of our lattice systems along the unit circle 



42 
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where k is the sequential index assigned to the zeros along the locus and is the angle 
subtended from the real-(/3J)^ axis to a given zero. The scaling of the zeroes density describes 
the renormalization group scale transformations of the lattice system, hence, 



A{(3, J, A, L) = L^A{(3L, JL, A) (20) 
where c is the scaling exponent in the vicinity of a second order fixed point. Therefore, 

^lim^A(J) = J^(l - J...) lnA(0) = 5/0+ ... (21) 

III. NUMERICAL RESULTS 

We have generated lattice ensembles of the spin-1/2 antiferromagnet on the square lattice 
at several different values of length L = {12, 16,20,24,32,40}, at the inverse temperature 
(3 = 100, at strong coupling (where J = 1) and at three different fixed values of the XXZ 
anisotropy parameter A = {1.01, 1.78, 2.0} using the continuous-time Quantum Monte Carlo 



method [37|. We have constructed the characteristic polynomial coefficients, defined in 
(12), for a polynomial ring whose singularities form the quotient using the numerical 
transfer matrix entries generated for the Markov process of each ensemble, following (18). 
We have applied standard rootfinding techniques to this polynomials to find its zeroes, and 
have constructed the simple difference relation for the zeroes density along the unit circle 
in the complex expansion parameter plane defined in (19). This zeroes density relation 
equivalently defines the Von Neumann entanglement entropy of the spin-1/2 antiferromagnet 
on the square lattice in (17), for the entanglement defined between the set of disjoint 
subsystems of the square lattice (that are analytic in the expansion parameter (3 J) and the 
rest of the ensemble. Each of these subsystems is therefore defined nonperturbatively by 
the mismatch between 0i ® 02 and 6i ® O2 in (4). This approach allows us to probe the 
entanglement entropy scaling of the primary chiral field associated with directly, even 
though a general state of the spin-1/2 antiferromagnet on the square lattice system is defined 
by two non-orthogonal polynomials [l^ . 

Whilst all of the zeroes we have evaluated in the complex plane automatically live on the 
surface traced out by 0i ® 02 (from the definition of the polynomial in (12)) the real test of 
this numerical method is whether the zeroes density (entanglement entropy) scales in any 
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meaningful way along this surface. As with the Ising model 4J], the zeroes polynomial in 
(12) defines the extension of the expansion parameter X{= (3 J), and this extension becomes 
the central extension when the lattice ensembles we have generated are both Lorentz and 
scale invariant. Hence, it is only if the zeroes demonstrate some symmetry by lying on a 
well defined locus in the complex plane that the extension of the system can be quantified 
via the scaling of the zeroes density. 

In Figure 1 the polynomial zeroes are plotted in the complex- (/3 J) ^ plane for three lattice 
ensembles generated at A = {1.01,1.78,2.0}. For the ensemble generated at A = 1.01 in 
Figure 1 all of the zeroes lie on the unit circle, hence the entanglement entropy scaling can 
be quantified via the zeroes density along this locus. However, whilst for A = 1.78 and 
A = 2.0 a subset of the zeroes lie on the unit circle, a further subset of the zeros lie on (or 
to the exterior of) a line at a fixed angle from the real-(/?J)^ axis. There is therefore a 
mismatch between 0i ® 02 and 6i ® 62 for these ensembles. The reason for this is that in 
(12) the polynomial is assumed to be analytic in A and J but it is only when the system 
is at criticality (Lorentz and scale invariant) that the scaling is critical in both variables. 
Hence, further from the Neel point scaling in A becomes subcritical and the critical point 
develops at = 0'. 

In Figure 2 we directly test the scaling hypothesis for 0i ® 02- We rescale the pseudo- 
critical value by subtracting from it the value of the zero closest the real-(/?J)^ axis 
00, and plot this difference (0^/ — 0o) against 1/L^. There should be no evidence of any 
meaningful scaling if the mismatch between 0i®02 and 61^62 is a purely numerical artifact. 
However, the result in Figure 2 indicates a clear linear confining flux between these two 
branch points, with a string tension of a = 68(4) for A = 1.78, and a = 78(7) for A = 2.0, 
in lattice units. The results in Figure 1 and Figure 2 strongly suggest that it is meaningful to 
model the general states of the spin- 1 / 2 antiferromagnet on the square lattice system either 



by two non-orthogonal polynomials 



I2I ]. or equivalently by a chiral wave theory coupled 



via a nonperturbative Chern-Simons term in the continuum limit ISj, since the polynomial 
zeroes we have measured do have a well defined extension, which links the UV to the IR via 
statistical transmutation jsljcl jl5|. 

On the basis of dimensional analysis 30| 3J], the general picture of entanglement entropy 
scaling in (2 + 1) dimensions is of a leading area law term with a prefactor of a^~'^ (where a is 
the lattice UV cutoff), which is therefore not universal, and a subleading universal area law 
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FIG. 1: The polynomial zeros of the spin- 1/2 antiferromagnet on the square lattice plotted in 
the complex- iiT^ plane (where K = /3J) for three different fixed values of the XXZ anisotropy 
parameter A = {1.01,1.78,2.0}. The lattice size of all three systems is = 32 x 32, and the 
inverse temperature is kept fixed at /? = 100. A subset of the zeros lie on a locus, corresponding 
to the unit circle in the complex-iT^ plane, for all three ensembles. 
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FIG. 2: Pseudo-critical scaling of the difference in the critical values of the zeroes, (/>o and (f)^i, as 
a function of the inverse lattice size squared, 

term propotional to (where ^ is the correlation length) of the form of (3). This general 
picture is therefore modified if the system is either gapped or gapless 5l| . In 3^ , conformally 
invariant systems at criticality have been found to satisfy this general picture. However, for 
the spinless gapped systems considered in 48], whilst the leading contribution is as above, 



the subleading behaviour differs in the quantum critical and non-critical regimes where it 
is either scales as in (3) or is a negative constant. In 49| 50|, it is argued that the general 
form of the scaling for systems with finite Fermi surfaces is for a leading term proportional 
to the product of the area between the subsystems and a logarithmic correction: a rescaling 
of the general picture on to the finite Fermi surface. Our expectation of the numerics from 
';hese results is therefore threefold: dependent on whether the lattice ensembles are gapped 
48 1 ■ gapless and conformally invariant 32] or gapless with finite Fermi surfaces j4^|50]. 

In Figure 3 the logarithm of the zeroes density in (19) is plotted as a function of the angle 
subtended along the unit circle in the complex- (/5 J) ^ plane, for an ensemble generated at 
A = 2.0. The scaling of the zeroes density can be fitted via two separate curves; where the 
density of the zeros sharply increases at the points ~ and cf) ~ 0', and also for the fiat 0- 
independent plateau. The scaling exponent values for fits to the second order scaling ansatz 
in (21) at ~ and ~ 0' are tabulated in Table 1, along with the value of the intercept 
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of the flat plateau, IhAtir. The fitted exponent values are different for ensembles generated 
at A = 1.78 (where c ~ 1) and A = 2.0 (where c ~ 1/9), whilst for A = 1.01 the measured 
scaling is purely first order. Although there is variation in the exponent values within the 
(jackknife) error estimates, it is important to note that the lattice units used to define 
in the analysis are not held fixed. The comparison of finite size effects through volumes 
is therefore best made through In Arm rather than L^. If the scaling exponent values are 
compared directly with comparable DMRG measurements the values show a similar level of 
consistency ClS db function of lattice volume [26l||27l]. 

For ensembles generated in the vicinity of the Neel point, at A = 1.01, the subsystem 
is analytically continued on to the boundary of the system 6i ® ^2- The fiat plateau in the 
zeroes density indicates that the entanglement entropy scales linearly in the boundary area 
(6*1 (g) 62) with a non-universal prefactor. There is no multiplicative logartihmic correction, 
and the scaling is therefore of the general form expected for conformally invariant critical 



points 



32l | 481 ] ■ Any subleading universal correction is lost in the analysis, because the 



analytic continuation is not defined beyond the system boundary. 

This scaling picture differs for the ensembles generated at A = 1.78 and A = 2.0, 
where the subsystem is analytically continued on to the boundary of the subsystem 
01 ® 02, and not on to the boundary of the system 61 ® ^2- For ensembles generated in the 
vicinity of the first order transition between the Ising phase and quantum phase above the 
Neel point, at A = 1.78 l^[l3], the flat plateau in the zeroes density indicates that the 
entanglement entropy scales linearly in the boundary area of the subsystem (0i ®02) with a 
non-universal prefactor and also logarithmically with a universal prefactor c. However, the 
01 ® 02 subsystem itself also scales linearly on to 61 ® O2 from the operator deflnition in (4) 
(this corresponds to the global mismatch between 0i ® 02 and 61 ^ 62). Hence the scaling 
follows the general form for gapless systems with a flnite Fermi surface in 491] 50|], of area 
law times a logarithmic correction. 

Even though the mismatch between 0i CS> 02 and 61 (g) 6*2 is a non-universal factor, by 
constructing exact operators for we are able to extract the central charge and Euler 
characteristic for the flnite Fermi surfaces of these ensembles. At A = 1.78 the ensembles 
lie in the vicinity of the flrst order transition between the Ising phase and quanturn phase 
above the Neel point, whereas at A = 2.0 the ensembles are in the Ising regime [l^[l3]- In 
both cases the central charge for the primary chiral fleld in is c = c = 1/2, but the Euler 
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FIG. 3: The logarithm of the zeroes density \nA{(j)) (= In X^ic/))) of the spin-1/2 antiferromagnet on 
the square lattice plotted as a function of the angle (j) subtended from the real axis of the complex- 
plane to a given zero along the unit circle. The lattice size of this ensemble is = 40 x 40, 
the inverse temperature is (3 = 100, and the XXZ anisotropy parameter is A = 2.0. 

characteristics differ as the later case the Z2 symmetry is broken. The Euler characteristics 
for the surfaces of the ensembles are defined by counting the number of vertices of the 
hypercubes in (8), subtracting the number of edges and adding the number of faces. At 
A = 1.78 the Euler characteristic is given by 12/3 — 8/2 + 6 = 6, whereas at A = 2.0 the 
Euler characteristics is given by. 8/3 — 8/2 + 2 = 2/3. From (3), the universal prefactors 
of the scaling of the boundary 0i ® 02 should therefore be c ~ 1 and c ~ 1/9, which are 
consistent with the values in Table 1. 

IV. SUMMARY 

We have presented a new procedure in this article, which uses the convergence properties 
of an exact expansion of the lattice partition function of the continuous-time Quantum Monte 
Carlo method to identify the critical scaling of the spin-1/2 antiferromagnet on the square 
lattice for ensembles we have generated at fixed XXZ anisotropy values of A = 1.01, A = 1.78 
and A = 2.0. This procedure is closely related to recent exact diagonalization studies of the 
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L 


A 


C 0^0 




InAr/R 


12 


1.78 


1.0332(0.1903) 


1.0781(0.0182) 


-2.3983(0.0090) 


16 


1.78 


1.0177(0.0226) 


1.0160(0.0256) 


-2.4309(0.0025) 


20 


1.78 


1.0097(0.0189) 


1.0074(0.0648) 


-2.4291(0.0021) 


24 


1.78 


1.0073(0.0250) 


0.9999(0.2222) 


-2.4236(0.0017) 


32 


1.78 


1.0101(0.0065) 


0.9997(0.0167) 


-2.4245(0.0012) 


40 


1.78 


1.0066(0.0480) 


1.0000(0.0119) 


-2.4052(0.0003) 


L 


A 






InAr/R 


12 


2.00 


0.1102(0.0952) 


0.1092(0.0526) 


-2.5679(0.0132) 


16 


2.00 


0.1502(0.0505) 


0.1547(0.0045) 


-2.3807(0.0022) 


20 


2.00 


0.1504(0.0471) 


0.1305(0.0034) 


-2.4685(0.0017) 


24 


2.00 


0.1401(0.0057) 


0.1076(0.0058) 


-2.4863(0.0010) 


32 


2.00 


0.1337(0.0067) 


0.1077(0.0065) 


-2.4690(0.0002) 


40 


2.00 


0.1226(0.0476) 


0.1172(0.0313) 


-3.5141(0.0038) 



TABLE I: Dependence of the scaling exponent c (determined from the fit of the logarithm of 
the zeros density to (j2ip ) on the lattice system length L and anisotropy A, at a fixed inverse 
temperature of /? = 100. In the final column we give the non-analytic first order contribution to 
scaling, InAr/ij, corresponding to the intercept of the flat plateau in Figure 2. 



scaling of the entang 



ement entropy of several quantum chain systems, and shows a similar 



level of accuracy 271] 26|] 34] • The entanglement entropy in our approach is found by calcu- 
lating the density of zeroes of a polynomial ring, and we have used this ring to represent the 
Chern-Simons phase difference between the two non-orthogonal polynomial representations 
that define the ground state of the spin-1/2 antiferromagnet on the square lattice in 12]. 
We have shown that the meromorphic convergence of the density of these zeroes arises from 
treating the replica method for entangelment entropy jsO] via a prescription for (^-function 
renormalization. We have compared our results with three different scaling pictures for 
the entangelement entropy of two dimensional fermionic quantum systems which are either 
gapped 48], gapless and conformally invariant [s^ or are gapless with finite Fermi surfaces 
491 ] SO]. Our approach allows us to identify the universality of the scaling associated with 



15 



the primary chiral field of tlie spin-1/2 antiferromagnet on the square lattice, even though 
the Fermi surfaces of our ensembles may be finite and the general form of the scaling of 
the entanglement entropy in this case comes with a non-universal prefactor 4^(5^. This 
scheme is therefore useful for directly investigating the entanglement entropy of spin impuri- 
ties in the underdoped cuprates that can be treated via the spin-1/2 antiferromagnet on the 
square lattice j?] 8]. Nothing in our construction forces the zeroes density we have obtained 
from numerics to scale in a meaningful way. However, in each of the ensembles we have 
generated, the polynomial zeroes show good evidence of an analytic continuation symmetry 
which validates the idea presented in 5| 61] 15|] that there is an intrinsic Chern-Simons term 
in the quantum regime of the spin-1/2 antiferromagnet on the square lattice that leads to 
the statistical transmutation of fermionic statistics from the UV into the IR. In the vicinity 
of the Neel point at A = 1.01, we have found that the scaling follows the form for gapless 
and conformally invariant critical points in 2d quanturn systems 32l ]. This supports the 
general idea presented in deconfined quantum criticality 21[ 22], that the Neel point can be 
smoothly deformed by quantum fluctuations to give a second order fixed point somewhere 
above A = 1. For our ensembles in the Ising phase at A = 1.78 and A = 2.0 we have we 
have found the scaling follows the form for gapless with finite Fermi surfaces in 2d quantum 
systems 49| 50|] . The Euler characteristics we have identified for the prefactor of the scaling 



associated with the primary chiral field of these ensembles are consistent with topology of 
the hypercubes we have used to construct the lattice operators in these two regions which 

nn 

are above and at the boundary for the Ising phase pjj p/7|] . 
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